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\  ABSTRACT 

^  Stationary  eigenmodes  are  derived  for  wave  propagation  in  a  medium  with 
quadratic  and  cubic  nonlinearities.  Aperiodic  algebraic  solitary  waves  are 
the  eigenmodes  with  a  continuous  spectrum  while  periodic  solitary  waves  are 
those  with  a  spectrum  comprising  an  infinite  number  of  harmonics.  Stationary 
eigenmodes  con^rising  the  fundamental  and  the  second  harmonic  only  have  also 
been  derived.  The  stability  of  the  aperiodic  solitary  waves  have  been  studied 
numerically  and  a  stabilization  technique  proposed.  The  stability  of 
stationary  eigenmodes  for  the  two-frequency  case  is  also  discussed."' 
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STATIONARY  EIGENMODES  AND  THEIR  STABILITY  DURING  WAVE  PROPAGATION 
IN  A  MEDIUM  WITH  QUADRATIC  AND  CUBIC  NONLINEARITIES  WITHOUT  DISPERSION 

P.  P.  Banerjee,  W.  Choe,  G.  Cao  and  W.  Hereman 

1 .  Introduction 

It  is  well  known  that  solitary  waves  can  exist  in  a  medium  with  a 
quadratic  nonlinearity  and  dispersion.  The  physical  basis  for  solitary  wave 
formation  in  such  a  system  is  the  critical  balance  between  the  steepening 
effect  of  nonlinearity  and  the  smoothening  effect  of  dispersion.  Since 
solitary  waves  do  not  change  their  shape  during  propagation,  they  can  be  termed 
stationary  traveling  wave  eigenmodes  (with  a  continuous  spectrum)  of  the  system 

[1]  .  [Unless  otherwise  stated,  we  will  use  the  term  "solitary  waves"  to  denote 
aperiodic  (e.g.  pulse-type)  solutions.]  If,  furthermore,  the  solitary  waves 
preserve  their  shape  upon  mutual  interaction,  they  may  be  called  "solitons" 

[2]  .  Periodic  solitary  waves  (e.g.  cnoidal  waves)  have  a  discrete  spectrum 
which,  however  contains  a  fundamental  frequency  and  an  infinite  number  of 
harmonics . 

Analogous  to  the  continuous  spectrum  case,  stationary  eigenmodes  may  also 
exist  for  the  case  of  a  discrete  number  of  propagating  frequencies.  As  shown 
in  [3],  stationary  eigenmodes  comprising  frequencies  (0^  and  2(0^  may  exist  in  a 
medium  with  quadratic  nonlinearity  and  dispersion;  and  perturbations  around  the 
steady  state  exhibit  periodic  or  "recurrent"  behavior  [4],  similar  to  the 
celebrated  Fermi-Pasta-Ulam  (FPU)  recurrence  [5] . 

In  this  paper,  we  look  at  the  possibility  of  stationary  eigenmodes  (with 
continuous  as  well  as  discrete  spectra)  in  a  medium  with  quadratic  and  cubic 
nonlinearities,  and  no  dispersion.  A  physical  picture  for  the  existence  of 
solitary  waves  may  be  portrayed  as  follows:  Consider,  for  instance,  a 
kinematic  wave  equation  of  the  form  [6] 

+  Cq  (1  +  ^ 

where  ^  represents  the  wavef unction,  c^  is  the  (linear)  phase  velocity 


and  P2'P3  denote  the  quadratic  and  cubic  nonlinearity  coefficients.  Consider, 
for  now,  the  case  where  ^  P3  where  y  at  t  ••  0  is  a  baseband 

(pulse-type)  signal  greater  than  zero.  Then  with  time,  the  leading  edge  of  the 
pulse  steepens  while  the  trailing  edge  smoothens  under  the  action  of  the 
quadratic  nonlinearity  alone,  while  the  reverse  occurs  under  the  effect  of  the 
cubic  nonlinearity  (see  Fig.  1) .  The  combined  effect  can  be  visualized  as  a 
balancing  process  whereby  the  signal  may  finally  evolve  into  a  shape  which 
remains  unchanged  during  propagation. 

The  organization  of  the  paper  is  as  follows:  In  Section  2,  we  discuss 
aperiodic  and  periodic  solitary  wave  solutions  of  a  model  equation  with 
quadratic  and  cubic  nonlinearities  only,  viz.,  the  nonlinear  Klein-Gordon 
equation  without  dispersion.  Results  for  one-dimensional  propagation  only  are 
presented;  for  extension  to  higher  dimensions,  the  reader  is  referred  to  [7] , 

In  Section  3,  numerical  simulations  on  interactions  of  aperiodic  solitary  wave 
solutions  are  presented.  Results  show  instability  upon  collision;  however, 
this  can  be  removed  by  incorporating  a  small  amount  of  "saturating"  fourth- 
order  nonlinearity  in  the  system.  In  Section  4,  we  study  the  two  frequency 
case,  viz.,  the  fundamental  and  its  second  harmonic,  by  starting  from  the 
)tinematic  wave  equation  (1)  and  establish  relations  between  the  spectral 
amplitudes  in  the  stationary  case.  A  stability  analysis  of  the  stationary 
eigenmode  in  the  discrete  frequency  case  is  also  performed;  these  results  are 
presented  in  Section  5.  Perturbations  around  the  steady  state  sometimes 
exhibit  periodic  behavior,  similar  to  FPU  recurrence.  For  the  solution  to  the 
initial  value  problem,  the  reader  is  referred  to  [8] . 

2 .  Aperiodic  and  Periodic  Solitary  Wave  Solutions  of  the  Nonlinear  Klein 
Gordon  Equation  without  Dispersion 

We  shall  take,  as  our  model  equation,  the  nonlinear  Klein-Gordon  (NKG) 


equation  without  dispersion  in  the  form 


(2) 


d^V  2  2  3 

72  ■  ^=0  72  ■  V  +  A3V  ;  V  =  V  (x,t), 

at  ax 

where  c^,  as  in  (1),  represents  the  linear  phase  velocity  and  where  A^,  A^ 
denote  the  quadratic  and  cubic  nonlinearity  coefficients.  In  our  quest  for 
solitary  wave  solutions  to  (2),  we  introduce  a  traveling  frame  of  reference  and 
renormalize  according  to 


c  = 


/■ 


"2  2 
c^  -  V 


2i 


1/2 


(x-vt),  A3  >  0, 


(3) 


to  get 


'f'  =  -  V  .  'J'  =  'i'(C) 


dC 


Now  multiplying  (4)  by  d'F/d^  and  integrating  w.r.t.  we  obtain 


dC 


12 


I  vjP  _  i  4,4  +  K  , 

3  2 


(4) 


(5) 


where  K  is  an  integration  constant.  Assuming  that  at  ^  =  0,  4' =  4'^  and 


d4'/dC  =  0,  we  get 

14  2  3 

^  i  '*'0  -  f  V  • 

In  order  to  reduce  (5)  to  a  tractable  integral,  we  put 

4'  =  4'  +  a 

and  set  the  constant  term  equal  to  zero;  this  gives 

.2 


(f)' 


(2a^  -  2a^)  4'  +  (2a  -  3a^)^  +  (|  -  2a)  4^  -  |  4'^ 


(6) 


(7) 


(8) 


with  either 


a  =  «!  - 

or  given  by  the  (real)  roots  of  the  cubic  equation 


(9) 


*  1%  - 1'  ♦  '■*'0  -  f  *  <V  f’ 


(10) 


In  what  follows,  we  will  discuss  the  two  cases  (a)  4/ 


0  3 


and 


(b)  “  3  tiy  or^e.  It  turns  out  that  case  (a)  yields  the  aperiodic 

(algebraic)  solution  while  case  (b)  gives  periodic  solutions. 


.1  Case  (a):  -  4/3 


The  possible  values  of  a  are  =  4/3  (from  (9))  and  =  0  (from  (10)). 
We  will  discuss  the  latter  possibility;  it  can  be  easily  shown  that  the  other 
choice  leads  to  identical  results. 

For  the  latter  choice,  K  =  0  and  the  ODE  in  (5)  can  be  readily  solved  to 
give  the  aperiodic  algebraic  solitary  wave  of  the  form  (see  Fig.  2) 


^'(C)  = 


1+2C  /9 


Note  tlie  similarity  of  (11)  with  the  solitary  wave  solution  of  the 
Benjamin-Ono  equation  [9] .  Note  that  the  amplitude  of  the  solitary  wave  is 
fixed  through  the  choice  of  the  nonlinear  coefficients  A^,  A^  (see  Eq. (3) ) . 
The  velocity  v(<Cq)  may  be  related  to  the  width  of  the  solitary  wave 
through  use  of  the  same  equation. 


2.2  Case  (b)  :  4'^  *  4/3 


The  possible  values  of  a  to  examine  are  now  given  by  (9)  and  the  roots  of 
the  cubic  equation  (10) .  We  state  here  a  lemma,  without  proof,  regarding  the 
roots  of  the  cubic  equation.  The  proof  of  this  is  detailed  in  [7]. 


Lemma  2 . 1  Of  the  three  roots  of  (10),  one  is  real  and  the  other  two  are 
complex  conjugates. 


Thus,  it  suffices  to  consider  the  two  real  roots  of  a(=a  )  for  each 

1,2 


value  of  'Fq.  Now,  from  (8),  upon  writing 


1 


(12) 


we  have 


where 


o  2  ^  3 
a  =  2a  -2a 


+  b<|i^  +  c4>  +  d  =  p(0), 

2  2 

;  b  =  2a-3a  ;  c  =  3*2a 


d 


1/2. 


(13) 

(14) 


The  solution  and  hence  '¥  and  'F,  can  thus  be  expressed  in  terms 

of  Jacobian  elliptic  functions  [10].  Hence  is  periodic  in  nature. 

Examining  (4),  it  is  obvious  that  if  d'P(0)/dC  =  0,  then  >  1(<1) 

corresponds  to  the  maximum  (minimum)  value  of  S'.  We  shall  call  the  maxim.um 

(minimum)  value  T. ,  ,  with  'F  ^  denoting  either  of  the  two. 

M(m)  ext 

Now,  it  may  be  readily  seen  that  V  and  a  both  satisfy  a  quartic 

GX  t 

equation  of  the  form  as  in  Eq.  (6),  with  replaced  by  ®rid  a 

respectively.  Since  we  have  shown  that  a  has  two  real  roots,  they  must  be  4^ 

and  S'  .  Thus  if  *F  ,  ,  is  the  given  initial  value,  y  will  follow  as  the 

m  M(m)  m(M) 

solution  to  (10)  .  Furthermore,  straightforward  arguments  may  be  advanced  to 

show  that  if  1  <  T  <  4/3,  then  0  <  S'  <1,  and  if  >  4/3,  then  T  <0. 

M  m  M  m 

In  order  to  explicitly  evaluate  the  solution  to  (13) ,  we  have  to 
calculate  the  roots  of  the  polynomial  p(<|>)  .  We  state,  again,  a  lemma  without 
proof  concerning  the  roots  of  the  cubic  polynomial.  The  proof  is  detailed 
in  [7]. 


Lemma  2 . 2  The  polynomial  p(4i)  in  (13)  has  only  one  real  root 
1 


P 


'F  -a 
ext 


(15) 


where  a  ,  if  'F 

M(m)  ext  m(M) 

The  above  lemma  exposes  four  different  cases  to  be  considered. 

Corresponding  to  a  given  initial  .condition  T.  »  4*  (='F.^,  .),  there  are  two 

0  ext  M(m) 

values  of  a;  a  >  1  and  o,  <  1 .  Suppose  the  initial  condition  is  S'  .  The  two 
1  ^  M 

values  of  a  are  'F,  and  4^  .  If,  next,  we  choose  the  initial  condition  as  the 

M  m 

'F  corresponding  to  the  4^  chosen  above,  we  will,  once  again,  get  two  values 
m  M 


of  (Xt  which  are  the  same  as  before.  Rigorous  relationships  between  these  four 
solutions  have  been  derived  in  [7];  in  this  writeup  we  will  only  present 
selected  results.  This  is  in  view  of  the  fact  that  the  second  order  ODE  as  in 
(4)  should  have  a  unique  solution  once  two  initial  conditions  (viz.,  the  value 
of  the  function  and  its  derivative  at  ^  0)  have  been  specified,  irrespective 

of  the  value  of  a. 


For  instance,  for  a  =  a,  “  ‘K,(>1)»  we  obtain  from  (13), 

i.  H 


where 


(16) 


1  I  3  *^1  2  ^1  *^1 

p  «{i)  =  -  p((t>)  =(|)  +  —  (J)+  —  (t)+  — 

^  ^1  a=a  ^1  ®i 

(using  (13)),  (18) 

where  a,,  b, ,  c  .  d,  are  defined  by  (14)  with  a  =  a,  =  ,  and 

1111  1  M 

denotes  an  integration  constant  to  be  determined  from  the  initial 
condition (s) .  After  some  algebra  and  using  (7),  (12)  and  (14)  it  follows 

that  [10] 

1  +  cn((2aJ-2aV^^  X  (C-^  )] 

'P  (C)  =a  +  - ^ ^ ^ -  (19) 


where 


=  lp[ 


and  where  the  parameter  cos  8^^  of  the  cn  function  is  given  by 


7  11  1/2 

cos  -  1  -  m^  =  -  +  -  Pj^”(Pj^)/IpJ^(Pj^)] 

The  period  of  the  cn  function,  and  hence  of  is  given  as  [10] 

^  ^  4K(l-mj^) 


(22) 


where 


K(4)  = 


n/2 


2  -1/2 

(1-H  sin  0)  d0. 


(23) 


Corresponding  to  initial  conditions  and  'f'  the  respective  constants 

M  in# 


in  (19)  may  be  calculated.  The  results  are,  using  (17) 


4' 

lM(m) 


(O 


+ 

M 


1  T  cn[(24'^^-2'F^^)^^^A.^C] 
(Pj+X^)  cn  (  (24^-24^) 


(24) 


It  is  easily  seen  that  there  is  a  half-period  shift  (=2K  (1-m  ))  between 

1  IM 

and  S',  ,  as  is  to  be  expected. 

Im 

As  an  example,  consider  the  case  where  4^^  =  2;  4'’  (0)  =  0.  From  (24),  the 
analytic  solution  reads 


«/  1  +  4.947  cn(1.530C) 

*  3.041  -  0.067  cn(1.530C) 


(25) 


with  1-m^  ==  0.571  and  A  =  5.015.  This  is  plotted  in  Fig.  3. 

In  passing,  we  remark  that  whether  the  periodic  solution  can  be 
decomposed  into  a  series  of  associated  aperiodic  solitary  waves,  as  can  be 
achieved  for  solutions  of  the  KdV  equation  [11],  is  still  under  investigation. 


3 .  Interaction  of  Algebraic  Solitary  Waves  and  their  Stability 

In  this  Section,  we  present  a  numerical  scheme  to  test  the  stability  of 
the  algebraic  solitary  waves  derived  in  (11)  upon  mutual  interaction  and 
propose  a  method  to  "stabilize"  them.  The  procedure  followed  is  similar  to  a 
recipe  from  nonlinear  optics  where  radially  symmetric  solutions  in  a  medium 
with  cubic  nonlinearity  are  sometimes  stabilized  by  incorporating  a  small 
amount  of  a  saturable  fifth-order  nonlinearity  in  the  system  [12] . 

The  numerical  scheme  essentially  enploys  an  explicit  finite-difference 


method  [13]  to  discretize  the  NKG  equation  without  dispersion  (Eq .  (2))  as: 


-  c^(AtMx)^  2y^)  -  (y^_^  -  2y^) 

+  (At)^  (A  (y^)^  +  A  (y^)^],  (26) 

23  -3  3 

where  i  and  j  refer  to  the  number  of  increments,  Ax  and  At,  in  space  and  time 
respectively.  However,  for  numerical  analysis,  it  is  required  to  show  that 
the  above  scheme  is  convergent.  Unfortunately,  the  proof  of  convergence  is 
rather  involved  and,  even  when  possible,  a  satisfactory  theory  can  only  be 
advanced  for  linear  PDEs .  We  follow  here  a  suggestion  by  Ablowitz  et  al  [14] 
that  a  certain  degree  of  stabilization  can  be  reached  in  the  numerical  scheme 
by  using  the  average  of  several  adjacent  space  increments  to  replace  y^  in  the 
nonlinear  term(3),  viz. 


(Vj)" — >  (J 


(27) 


and  by  using  Ax  = 

Fig.  4  shows  the  three  regions  we  employ  for  our  computer  simulations.  I 
and  II  are  linear  nondispersive  while  III  represents  the  quadratically  and 
cubically  nonlinear  region.  Our  numerical  experiment  involves  malcing  two 
(identical)  pulses  traveling  to  the  right  and  left  in  regions  I  and  II 
respectively.  These  pulses  eventually  enter  region  III  where  they  interact 
with  each  other.  Outflow  (matched)  boundary  conditions  are  simulated  at  x  =  0 
and  X  =  8  (see  Fig.  4)  to  prevent  any  reflections  from  the  extremities. 
Interface  boundary  conditions  at  x  =  2  and  x  =  6  are  set  such  that  the 
wavef unction  and  its  spatial  derivative  are  continuous  across  each  interface. 

The  initial  conditions  are  set  from  a  (knowledge  of  the  single  aperiodic 
solitary  wave  solution.  Note  that  from  (11),  the  denormalized  algebraic 
solitary  wave  may  be  written,  using  (3),  as 


1  + 


4A_ 

_ ^ 

3A, 


2A, 


9A. 


(x-vt) 


7 


) 


y(x, t) 


(28) 


As  one  initial  condition  we  thus  set 


y(x, 0)  =  2  X  (-4/3) (A2/A2) 


1  +  ^4 

2  2 

9A3(Co-v  ) 


1  +  ^4 

9A3(c^-v2) 


(29) 


where  x^,  x^  lie  in  regions  I  and  II  respectively.  The  differential  initial 
condition  viz.,  y^(x,0),  is  set  equal  to  zero.  This  facilitates  propagation  of 
each  pulse  in  both  positive  and  negative  directions,  starting  from  each  linear 
region.  At  t  =  0,  the  right  and  left  propagating  pulses  overlap,  accounting 


for  double  the  amplitude  2 


4A, 


3A. 


After  some  time,  each  initial  pulse 


breaks  up  into  two  with  amplitudes  equal  to 


one  traveling  toward 


region  III,  the  other  toward  an  extreme  matched  boundary  (at  x  =  0  or 

X  =  8)  . 

2  2 

We  choose  v  =  1.0,  c^  =  1.0050  [c^  -  v  =  0.01],  A^  =  -1,  A3  =  1  and 
^1  =  1,  X2  =  6  for  our  numerical  experiment.  The  results  are  plotted  in  the 
sequence  of  Figs.  5 (a) -(g),  from  which  we  may  easily  see  that  the  initial 
waveforms  distort  severely  due  to  interaction  in  the  nonlinear  region.  No 
appreciable  distortion  is  observed  at  the  linear-nonlinear  interface. 

In  what  follows,  we  show  that  the  observed  distortion  may  be  removed  by 
incorporating  a  higher  order  saturating  nonlinearity  in  the  system.  The 
equation  describing  the  modified  system  can  be  expressed  as 


dt^ 


afy 

ax^ 


A2'*' 


A3V 


A^y 


<  0, 


0,  A^  <  0.  (30) 


Eq.  (30)  ia  then  programmed  with  the  same  values  of  A^/  A^,  v  amd  as 

before,  and  with  A,  =  -0.12.  This  value  of  A.  was  chosen  from  a  series  of 
4  4 

trial  values  incorporated  to  minimize  distortion  upon  interaction.  The 
results  are  plotted  in  the  sequence  of  Figs.  6 (a) -(g),  showing  that  with  the 
help  of  stabilization  by  the  saturating  fourth-order  term,  the  algebraic 
solitary  waves,  again,  preserve  their  shape  upon  interaction.  Numerical 
experiments  with  two  solitary  waves  moving  in  the  same  direction  (but  with 
different  velocities,  and  hence,  widths)  and  with  sequences  of  algebraic 
solitary  waves  are  currently  under  way. 

4.  Steady  States  in  the  Discrete  Case:  Fundamental  and  Second  Harmonic 


In  the  discussion  above,  aperiodic  and  periodic  solitary  wave  solutions 
were  derived  for  a  medium  with  quadratic  and  cubic  nonlinearities,  using  the 
NKG  equation  without  dispersion  as  a  model.  The  aperiodic  solitary  waves  can 
be  regarded  as  stationary  traveling  wave  eigenmodes  but  having  a  continuous 
frequency  spectrum.  The  periodic  solitary  waves,  however,  can  be  thought  of  as 
a  collection  of  discrete  frequencies  (viz,  the  fundamental  and  its  infinite 
set  of  harmonics)  that  form  a  stationary  eigenmode.  In  this  section,  we  will 
show  that  it  is  possible  to  establish  stationary  eigenmodes  even  for  two 
frequencies,  viz.,  the  fundamental  and  its  second  harmonic,  in  a  medium  with 
qfuadratic  and  cubic  nonlinearities.  We  will  ta)^e,  as  our  starting  point,  the 
nonlinear  l^inematic  wave  equation  as  in  (1),  which  though  different  from  the 
NKG  equation  (2),  nevertheless  provides  a  good  m.odel  for  such  a  medium. 

To  study  the  spatial  evolution  of  the  spectral  amplitudes  at  the 


fundamental  (0)^)  and  the  second  harmonic  (2(0^),  we  substitute 


1  T— > 

\|f(x,t)  “  2 


jn  (w^t-jc^x) 


in  (1)  with  G, 


••  G  to  ensure  that  the  wavefunction  is 
-n  n 


0  and  G 


real.  The  ratio  to^/kQ,  where  k^  represents  the  propagation  constant  of  the 
fundamental,  is  taken  to  be  equal  to  the  linear  phase  velocity  c^  in  (1) .  In 
general,  the  propagation  constant  for  the  nth  harmonic  may  differ  from 
nk^;  however,  if  so,  it  should  show  up  in  the  expression  for  (x) .  Using  the 
"slowly-varying"  assumption 


- —  «  nk„G  , 

dx  On 


together  with  the  normalizations 


A  ‘^*^2  A 

ue  =  G^/R,  ve  =  G^/R  , 


R  =  (  G^<0)  +  G^{0)  ) 


E  =  r  6_Rk.x, 

■9  2  ^2  0 


and  the  definitions 


A  M 
"  2P2  ' 


e  =  02  -  20^, 


we  finally  arrive  at  a  set  of  three  coupled  equations  linking  the  (real) 

spectral  amplitudes  u,v  and  the  relative  phase  difference  0. 

du  .  „  (35 

—  *  uv  sin  0 . 
dv 


uv  .  n 

-r  »  -u  sin  e. 


(35b) 


(2v  -  — )  cos  0  +  2Q(v  -  u  )  , 

V 


(35c) 


with  the  relation 
2  2 

u  +  V  ”1  (36) 

expressing  the  conservation  of  energy. 

Eqns .  (35)  may  be  compared  with  the  evolution  equations  derived  by 
Armstrong  et  al  [14]  in  connection  with  harmonic  generation  in  optics  in  a 
(quadratically)  nonlinear,  dispersive  dielectric.  While  Eqns.  (35a, b)  are 
identical  to  Armstrong's,  in  (35c)  only  the  first  term  on  the  RHS  is  identical, 
while  the  second  term  is  to  be  compared  with  As  in  [14]  which  represents  the 
phase  mismatch  between  the  fundamental  and  the  second  harmonic  in  a  dispersive 
medium. 

The  "relative  phase-loc)ced"  stationary  eigenmode  may  now  be  defined  as 
follows : 

du  dv  d0  -  (37) 

55  ■  5?  -  5f  ■  “ 

From  (35c),  this  means  (using  (33)  and  (34),  (36)) 


Is 


2  +  alG 


2sl 


1  +  aiG 


1/2 


2s 


A 

;  a  = 


P3  ■ 


23 


^2  cos  0 


(30) 


where  the  subscript  s  has  been  used  to  denote  the  steady  state. 

Observe  from  (35a, b)  that  possible  values  of  0  in  this  case  are  0  and  7t  (cos 

3 

0  =  ±  1) .  The  variation  of  the  second  harmonic  amplitude  with  the  fundamental 
s 

amplitude  for  the  stationary  eigenmode  is  plotted  in  Fig.  7.  Note  that  when  a 

<  0,  |G  I  is  a  multi-valued  function  of  |G  I .  The  explicit  relations  for  0 

^3  X  S  X 

and  obtained  from  the  respective  explicit  differential 

equations  (not  written  here)  which  occur  as  an  intermediate  step  in  the 

derivation  of  (35c)  .  It  may  be  readily  checited  that  <[1^^  and  vary  linearly 


with  4-  Thus,  the  propagation  constants  of  the  fundamental  and  the  second 


harmonic  are  indeed  modified  in  the  steady  state.  However,  the  new  (constant) 
phase  velocity  is  different  from  c^ . 

In  the  following  section,  we  will  perform  a  stability  analysis  of  the 
stationary  eigenmode (s)  to  find  which  branches  of  the  steady  state  plots  (see 
Fig.  7)  are  stable  and  which  are  not.  Whenever  possible,  the  period  of  the 
perturbation  will  also  be  calculated. 

For  a  description  of  the  solution  to  the  initial  value  problem  starting 
from  Eqns.OS),  the  reader  referred  to  [8]. 

5 .  Stability  Analysis  of  the  Stationary  Eiqenmode  (Fundamental  and  Second 
Har  monic) 

We  first  rewrite  (35c)  as 

— =  (2v-—  )  (l-3in^9)  ±  2Q(v^-u^)  v/l-sin^G  (39) 

dg  V 

and  perturb  the  system  (35a, b),  (39)  as: 

u  =  u  +  Au  , 
s 

V  =  V  +  Av  , 
s 

sin  0  «  (sin  0)  +  A(sin  0) .  (40) 

s 


Up  to  first  order  in  the  perturbations,  this  gives  the  following  linearized  set: 


d  (Au) 


u  V  A (sin  0) , 
s  s 


d(Av)  2.  .  . 

— —  =  -u  A(szn  0), 

dg  s 


d[A(sin  0)  ]  ....  -2  .  , 

— ‘ - -r - =  (3  +  V  ±  8Qv  )  Av. 

dg  s  s 


(41) 


Now  assuming  the  explicit  forms  of  the  perturbations  as  proportional  to 
e  the  permissible  values  of  \  (for  nontrivial  solutions)  are: 


0, 

»±jvu^(3  +  v  8Qv  )  . 

■'  a  s  a 


(42) 


The  condition  for  stability  requires  that  the  expression  under  the  radical  sign 
should  be  positive;  in  denormalized  form  this  translates  to 


4R^-3|G  1^ 

- - - r  for  a  ^  0. 

4a{R  -|G  I  ) 
is 


(43) 


Consider  the  case  a  >  0.  It  may  be  easily  chec)ced  from  (43)  that  any  value  of 

IG,  I  and  |G_  i  (as  given  in  Fig.  7(a))  will  be  stable.  For  a  <  0,  the  stable 
Is  2s 

and  unstable  regions  are  drawn  in  Fig.  8.  When  this  figure  is  juxtaposed  on 
Fig.  7(b),  it  is  clearly  seen  that  the  upper  "branch"  in  Fig.  7(b)  is  unstable, 
while  the  lower  "branch"  is  stable.  We  conclude,  therefore,  that  the 
stationary  eigenmodes  in  a  medium  with  quadratic  and  cubic  nonlinearities  can 
have  two  stable  states  in  the  two-frequency  (fundamental  and  second  harmonic) 
case . 


6 ,  Conclusion 

Stationary  traveling  wave  eigenmodes,  with  continuous  as  well  as  discrete 
spectra,  have  been  derived  for  a  medium  with  quadratic  and  cubic 
nonlinearities,  and  no  dispersion.  The  physical  basis  for  the  existence  of  a 
balance  between  the  different  nonlinear  effects  has  been  established. 

Aperiodic  (algebraic)  solitary  waves,  which  occur  as  solutions  to  the  NKG 
equation  without  dispersion,  are  an  illustration  of  stationary  eigenmodes  with 
continuous  spectra.  These  are,  however,  unstable  upon  mutual  interaction,  but 
may  be  "stabilized"  by  adding  the  right  amount  of  a  higher  order  saturating 
nonlinearity  in  the  system.  Periodic  solitary  waves  in  such  a  system  may  be 
expressed  in  terms  of  ratios  of  'cn'  functions,  and  may  be  loo)ced  upon  as 
stationary  eigenmodes  with  discrete,  albeit  infinitely  many,  frequency 
components.  For  yet  another  nondispersive  system  which  is  quadratically  and 
cubically  nonlinear,  but  modeled  by  a  different  equation,  we  have  established 
stationary  eigenmodes  comprising  the  fundamental  and  the  second  harmonic  and 


examined  their  stability.  Not  reported  in  this  paper,  but  presented  in  (8)  is 


the  fact  that  in  the  solution  to  the  initial  value  problem  in  this  case. 


periodic  exchange  of  energy  between  the  fundamental  and  the  second  harmonic  is 
observed,  similar  to  FPU  recurrence. 
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